Electronic structure and magnetism in slightly doped SrB 6 . 



O 

o 



oo 



o 

CZ3 



T. Jarlborg 

DPMC, University of Geneva, 24 Quai Ernest- Ansermet, CH-1211 Geneve 4, Switzerland 

(February 1, 2008) 

Spin-polarized band calculations for supercells of SrBg, where a La-, In- or Al- impurity or a vacancy 
is replacing one Sr, are performed within the local spin density approximation. Moderately large 
cells with 8 formula units (56 atoms) are studied for all dopings and large ones with 27 formula units 
(189 atoms) for the case of La-doping. The undoped system has a vanishing density-of-states (DOS) 
at the Fermi energy (Ef), while the additional La-d band makes Ef to enter the bands above the 
gap. An Al (or In) impurity has the opposite effect, with a rigid-band like shift of Ef to below the 
gap. In the former case, the addition of a d-electron makes a local, impurity-like modification of the 
electronic structure close to the La atom. As has been shown in a previous publication [Phys. Rev. 
Lett. 85, 186, (2000)], it can lead to a weak ferromagnetic state. This result is shown to persist 
at even lower doping, by the 189-atom cell calculations, with a moment of the order 0.1 \ib per 
La impurity. An addition of a p-electron, by an In-impurity, makes no similar effect. The DOS at 
Ef and the Stoner factors are in all cases low, and would not suggest ferromagnetism alone. The 
magnetic state can be understood from a charge transfer and an additional gain in potential energy, 
as a spin-splitting is imposed. A DOS with a localized impurity band is essential for this effect. 
Different models with modified localization of the La-band are made in order to show the correlation 
between charge transfer and the size of the magnetic moment. 
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I. INTRODUCTION. 



The structure of hcxaborides, XBg, is simple cubic, 
with one large atom (X = Eu, Ca, Sr or Ba) in the cor- 
ner, and a Bg octahedron in the center. It can be viewed 
as the CsCl-structure, where the Bg-octahedron would 
be a pseudo-atom in the center and the X-atom in the 
corner is the second atom. The internal parameter, z, is 
a measure of the size of the octahedron. It is the short- 
est distance from one of the B-atoms at the edge of an 
octahedron to the cube surface. Its value is about 0.2 of 
the lattice constant, but it varies from system to system. 
The hexaborides have been studied mostly because their 
intriguing properties showing signatures of both metal- 
lic and semi-conducting behaviorErtl. Recent experimen- 
tal works on La-doped hexaborides like La x Can_ x )B6 
show unexpectedly a very weak ferromagnetic state that 
persists up Xo large temperatures, and for very low x, 
about 0.005u. Independent measurements confirm these 
findings in part, but they also show that the results can 
span from weak ferromagnetism to strong paramagnetism 
and that Al impurities may conceal some of the intrinsic 
properties^. 

Ferromagnetism in the La-doped hexaborides XBg, 
X=Ca,Sr or Ba, is surprising because the density-of- 
states (DOS) at the Fermi energy, N(Ep), and the Stoner 
factor, S, are expected to be small in this family of com- 
pounds. Band calculationsQ in the local density approx- 
imation (LDA)Ej on SrBg and CaBg confirm that the 
DOS is low, with Ep falling in a valley of the DOS, and 
spin-polarized calculations for EuBg show a very weak 
spill-over of the spin density from the localized Eu-f spin 
to the rest of the valence band. The exact behavior 
of the band overlap at the X-i»oint has been monitored 
in several band calculations!!]. It has been found that 



this overlap, which determines if the undoped system is 
metallic or semi-conducting, depends on the structural 
parameter z. In fact, the overlap is more sensitive to 
z-variations than to variations of the lattice constant!]. 
A recent band calculationlS determined the fine details 
of the bands near Ep and the Fermi surface (FS), and 
concluded that they are in good agreement with de-Haas- 
van-Alphen measurementsE3. An experimental determi- 
nation of the._b.and dispersion the near the X-point, by 
photoemissonEa, is in fair agreement with the theoretical 
predictions for the existence of a X-pocket in the bulk, 
but a different interpretation is made for bands in the 
surface region. 

Electronic structure calculations based on the local 
(spin) density approximation, either by spin-polarized 
calculations or by calculations of the Stoner factor S, 
are in [general able to detect if a system is ferromagnetic 
or nolo O Because the Stoner factor_ior La^Sr^^Bg 
is far below the limit for magnetisnxH, it is tempting 
to suggest that the observed magnetism is not a band 
magnetism, but something else. Electron-hole interac- 
tion for an asymmetric position of Ef relative to the 
gap has been proposed to lead to excitonic magnetismlill. 
Spontaneous polarization of a dilute electron-gas has also 
been suggestedtrEa. However, for Ce it is found that the 
Stoner factor underestimates the tendency towards mag- 
netism, because the stabilization energy for a magnetic 
state comes also from the Coulomb energy in addition 
to the exchange energyo. The effects can be larger in 
a compound as long as the different sub-bands and the 
ionicity of the system show favorable conditions. In a 
recent work we showed that LaSr7B48 is likely to be such 
a system!!.. Here we extend these calculations to larger 
unit cells, to different types of dopants, and give a more 
complete discussion of the results. 
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II. THEORY. 

A. Band method. 

Electron energies are calculated by the linear muffin- 
tin orbital (LMTO) methodic using the LDA. The in- 
ternal parameter z, is in these calculations 0.215 of the 
lattice constant. The lattice constant is fixed at 7.93 a.u. 
in each calculation. Additional 'empty' spheres are in- 
serted at the cube edges, between the Sr-atoms, in order 
to account for non-sphericity of the potential in the open 
part of the structure. Calculated results are presented 
for two unit cells: One supercell of 8 elementary cubes (2 
cube lengths in x,y and z-direction) , containing totally 
80 sites, 56 atoms and 24 empty spheres. A second set of 
calculations was made using a larger super cell contain- 
ing 27 elementary cells (3 cube lengths in each direction) , 
with totally 270 sites (189 atoms plus 81 empty spheres). 
These two cells allow to study doping concentrations x 
of 0.125 and 0.037, respectively. Nonequivalent sites are 
treated independently. One basis set for the 80-site cell 
unit cell includes s,p,d for the atoms, and s,p for empty 
spheres. The results from these calculations can be found 
in refJld. A reduced basis set for the large 270-site calcu- 
lations includes s and p states for the B-atoms and only 
s-states for the empty sites in the direct basis. We include 
the B-d and empty-p "tail" states in the reduced basis 
calculations. This means that the d-states for B atoms, 
for example, are constructed from tail states coming from 
neighboring sites. This gives a good description of tail 
states like B-d and p-states on the empty sites, close to 
what was obtained from the large basis, without having 
to increase the size of the eigenvalue problem. Intrin- 
sic "on-site" states, like the high energy La-f state or B-d 
states, can not be constructed from tail states in this way, 
but these states are unoccupied and not important in this 
study. Test calculations for the smaller of the cells using 
this smaller basis showed very similar results as the ones 
from using the large basis, when the linearization ener- 
gies (E^) of the La-band are increased upwards about 0.2 
Ry. 

In the 270-site calculations, where we employed the 
same basis (with tail states) and the same E/s, which 
in the small basis 80-site calculation reproduced well the 
results of the older large basis calculation. In this way 
we believe that the results for the large 270-site cell are 
realistic despite the use of a minimal basis set. No lattice 
relaxation around the impurities is considered. The num- 
ber of k-points varies from 10 to 20 irreducible k-points 
for the unit cell of 80 sites, and from 4 to 10 for the 270 
site cell. 



B. Stoner models. 

The original Stoner modclEl can be applied to- 
gether with ab-initio paramagnetic band results in or- 



der to studj«_ .the spin susceptibility or the onset of 
magnetismli3ii3. An applied exchange splitting, £, leads 
to a transfer from minority spins to majority spins, and 
a loss in kinetic energy AK = N£ 2 , where N, the DOS 
at Ep, is assumed to be constant. The same spin trans- 
fer leads to a gain in exchange energy AE X = N 2 I S ^ 2 , 
where I s is the exchange integral. If a partial DOS is 
not constant near Ep, there will be a possibility of charge 
transfers, Aq^ and AQ, as function of £ and temperature, 
which make additional contributions to the total energy. 
It is possible to identify one sub-band, the La-impurity 
(mainly d) band, which has a large derivative of its DOS, 
N' La (E F ). The total DOS at E F is N = N La + N v , 
where N v is the rest of the DOS. By neglecting all other 
derivatives one can make a model of the charge transfers 
as function of £ and ksT . By taking into account the 
energy associated with the charge transfer, Uq, one ob- 
tains a modified Stoner criterion. This is made in detail 
previouslyEj and the new Stoner criterion is: 

Su = NI S + U N' La N v /N 2 > 1 (1) 

This modified Stoner criterion is useful for an under- 
standing of why a system can be pushed towards mag- 
netism by Coulomb energies. But, as will be discussed in 
section 3.2, the effects are too strong to use eq. [TJ for a 
quantitative prediction of a magnetic state. 

III. RESULTS. 

A. Bandstructures. 

The bandstcucture for undoped SrBg agree well with 
the one in ref.0. The Fermi energy falls in a low DOS 
region between the B-p and Sr-d band regions, where a 
small gap can appear. The B octahedron behaves partly 
as one identity, as a large pseudo atom. This is seen in 
a narrow band about 0.9 Ry below Ep. There is only 
one such band per basic unit cell, although in reality it 
is of B s- and p-character. This band is separate from 
the other B s- and p-bands. Its presence could be viewed 
as one single s-band from each Bg pseudo atom within a 
CsCl-structure, where Sr is one atom and Bg the other 
one. However, this is far below Ep and this way of look- 
ing at the B-octahedron is probably not useful for the 
bands near E^ or for the physical properties of doped 
hexaborides. 

The effect of doping, where one Sr is exchanged with 
another atom, show some unexpected result depending 
on the dopant. Fig |l| shows the DOS for undoped, for 
La-doped and In-doped SrBg, all for a;=0.125 (small 80- 
site cell, small basis). La has a d-band coming in just 
above Ep. In this energy range there are very few B and 
Sr bands for the La-d band to hybridize with. The La-d 
band can hybridize between La neighbors if the La con- 
centration is large, but if there are only few La atoms, 
the La-d band will localize and becomes narrow. Indium 
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has one additional valence electron compared to Sr, ex- 
actly as La has. But the additional electron in In is a 
p-electron, and in contrast to the d-band for La, the In-p 
band hybridizes well with the B-bands below ~Ep. The 
new In-band below Ep can accommodate two electrons, 
but In has only one additional electron compared to Sr. 
The result is that Ep will move downwards compared 
to the rest of the bands. Thus, as seen in fig. [y, the 
effects of In and La doping are opposite. Furthermore, 
with La there is a large increase of the local La-DOS, 
while the DOS on the In-site and a Sr-site are very sim- 
ilar. The resulting DOS from In doping is very simi- 
lar to the one from Al (or even B) doping, at least as 
long as lattice relaxation can be ignored. It has been 
reported that Al can penetrate the hexaboride material 
when the aluminum-flux-growth technique is used to pre- 
pare single crystalsO. The calculated results for single In 
or Al impurities shown here, indicate that the effects of 
Al contamination should be very different from those of 
La impurities. The effects could even appear quite unde- 
tected compared to undoped hexaboride, since the main 
result is a rigid-band like downward shift of Ep on a uni- 
form DOS. Theories of excitonic magnetism are based .qji 
rigid-band assumptions of the undoped band structureO. 
Such theories should therefore apply to Al (or In) dop- 
ing, but not easily for La-doping, since the latter is not 
rigid-band like. However, according to our results, the 
effect of the p-electron doping in this material is oddly 
transformed into hole doping. The charge distribution is 
such that about one half of the additional electron from 
In stays on the site, a bit more than one half is added 
close to the impurity. At the most distant sites there is 
slightly less electronic charge compared to the values for 
the undoped case, which presumably is an effect of the 
shift of Ep . 

The paramagnetic (non-spinpolarized) DOS from the 
calculation for LaSryB48 is shown in fig. [j]. The La band 
forms a separate band near Ep below the rest of the 
Sr-d band, making the real space distribution to be lo- 
calized near the La impurity from this band. The total 
LaSr 7 B 48 DOS at E F ,N(Ep), about 48 states/cell/Ry, 
is concentrated on the single La site, with NL a (Ep) ~ 
19 states/cell/Ry, cf Table 1. On the remaining 55 
sites, each Sr has on average 1.0 and each B roughly 
0.4 states/cell/Ry, respectively. By comparison, a rigid- 
band model obtained by adding one electron to the band 
structure without the La impurity gives very different 
values: Of the total DOS, 32 states/cell/Ry, each Sr has 
1.4 and each B 0.25 states/Ry. Clearly, La behaves like a 
localized impurity site with a very large local Ni, a (Ep), 
much larger than on Sr. The DOS values are included in 
Table 1. 

The additional valence electron remains located on the 
La. By comparing the charges in the undoped and doped 
supercell one finds that the remaining Sr7i?4g-spheres re- 
ceive no additional charge, cf. Table 2. The difference in 
charge between the two systems is never larger than 0.01 
electron per Sr or B site. The derivative of the DOS be- 



comes very large, about 1400 states/cell/Ry 2 , from which 
700 comes from the La site, compared to about 30 for 
each Sr site when the rigid-band model is applied to the 
undoped bandstructure. This asymmetric and uneven 
distribution of the DOS near Ep is a condition for real 
space charge transfers as T increases and as an exchange 
splitting is imposed. The transfer goes in both cases to 
the La site from the rest of the system. 

Next, for the calculations with the large 270-site 
cell only the case with a La-impurity was considered, 
LaSr26Bi62- The total and partial La-DOS near E^ are 
shown in fig^[ together with the corresponding DOS from 
the 80-site calculation. Both results are calculated with 
the same basis. The DOS functions from the two calcu- 
lations are divided by 27 and 8, respectively. This makes 
the units of the DOS states per formula unit (f.u. of 
SrBg), and the two cases can be compared on the same 
diagram. The number of k-points, 20 for the small and 
10 for the large cell, reflects a larger density of points in 
the reciprocal space for the 270-site calculation. It is seen 
that the La-band is more separated from other bands and 
more localized in the large cell. The value of the total N 
at Ep (6 states per Ry and f.u.) is quite similar in the two 
cases, but the local DOS on the La-site is much larger in 
the largest cell. The 270-site cell has a large number of 
sites far from the impurity with low DOS, almost as for 
the undoped material. The number of such sites is com- 
paratively larger than for the 80-site cell. But the very 
large DOS at the impurity site compensates for the low 
DOS elsewhere, and the average DOS per f.u. becomes 
similar for the two cells. The local La-DOS, which in the 
small cell is about 19, is 55 in the large cell, cf Table 3. 

A comparison of the DOS near E^ with the ones from 
the small cell, reveals more rapid variations of the DOS 
structures. The La-band is narrower, with the total DOS 
going from zero to very large values within a much nar- 
rower energy interval than in the case of the 80-site cell. 
The effect is larger DOS derivatives, mainly steming from 
La, but also larger sensitivity of the charge transfers and 
their T-dependencies. The amount of charge transfer 
as function of ksT (and £) is typically 3-5 times larger 
within the large cell, but the non-linearities are large and 
a small shift of E^ can make a large difference. It is 
therefore not meaningful to quote values of N' in this 
case. At the most distant (from the La) B sites, the 
local B-DOS is small, but still about 50 percent larger 
than what is the case for the undoped DOS, when a rigid 
band shift has been done to account for one additional 
electron. For the distant Sr-sites, the local DOS tend 
to be close to the value of the undoped DOS with the 
proper rigid-band shift of Ep, about 0.7 per Ry. The 
first 1-2 layers of Sr and B surrounding the La impurity, 
have approximately doubled their local DOS compared 
to the undoped values. 

Despite the large number of atoms separating the La 
in the largest cell (188 atoms, 3 types of Sr and 9 types of 
B), and despite the fact that the additional La electron 
stays localized near the impurity, there is still an influ- 
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ence on the electronic structure far from the impurity. 
Their DOS values are larger than the ones of the undoped 
DOS, which are obtained after a rigid-band shift of Ep 
to account for an additional electron. An even lower im- 
purity concentration seems to be required for finding the 
electronic structure far from the impurity, almost "insu- 
lating" like for the undoped case. In such a case, the 
interaction between local moments would effectively be 
cut-off. 



B. Magnetism. 

We first present the results using the 80-site cells, 
mostly using the small basis set as explained above. Spin- 
polarized calculations for undoped SrBg do not give sta- 
ble moments of significant size. An applied field is used 
to get a moment of about 0.5 \xb per cell, but this mo- 
ment decreases steadily during the iterations when the 
field is removed. The same situation holds for InSryE^g. 
A number of calculations were also attempted where one 
Sr was replaced with an 'empty' site. This was motivated 
by the reports that a kind of " self-doping" , by Sr vacan- 
cies might be present because some-,seemingly undoped 
samples show the weak magnetisms The unpolarized 
calculations for the cases with a Sr-vacancy showed the 
effect of charge transfer as function of T, so that as T 
increased there was a larger charge on the remaining Sr 
sites. However, the moments always decreased (to be- 
low 0.005 /is) in the spin-polarized calculations, also for 
the case where a slight B-relaxation towards the missing 
Sr-site was attempted. Thus, in the following we discuss 
only La-doping, since undoped, In- (or A1-) doping, and 
Sr-vacancy "self-doping" all failed to show magnetism. 

The 80-site cell calculation of LaSr7B4g using 20 k- 
points and the large basis, finds a small moment of 0.10 
/^s/cell, which is about 10 mRy/superceLLlower in to- 
tal energy than the paramagnetic solutior£3. If fewer k- 
points are used, one finds that the (paramagnetic) DOS 
and its energy derivative near Ep are larger, cf Table 1. 
One case with the small basis is also shown. The mo- 
ments are not quite as large as with the large basis, but 
qualitatively the results are the same. 

A number of different calculations with the 80-site cell 
and the small basis are done where the linearization en- 
ergies are varied in order to obtain different values for 
N(Ef) and its derivative. All these different cases are 
not presented in detail, but the trend is that when the 
local La-DOS and its derivative at are large, then 
the moments will be large (larger than about 0.1 [Ib per 
cell) in the spin-polarized calculations. Thermal broad- 
ening will usually oppose a large moment. But there are 
a few exceptions, when the energy derivative of the DOS 
is varying around E^?, showing that the sharpness of the 
DOS is not always close to Ep. Temperature dependent 
charge transfers also show deviations from a simple T 2 
behavior. The moment is always largest on La, with more 



than 50 percent of the total moment on and near the La 
impurity. The moment on B-sites far from the impurity 
is negligible. This fact and the uneven distribution of 
local DOS show that a rigid-band picture of dilute La- 
doped hexaborides does not apply. The bands near Ep 
are modified by the spatially non-uniform perturbation 
of the potential around the impurity site. 

The calculated Stoner product NI S is 0.210, N' La ~ 
730 per Ry 2 - cell, N - 48 and N v ~ 29 (per Ry • cell), 
to give w 440£ 2 and AQ = 720(k B T) 2 , with £ and 
ksT in Ry. With Uq calculated from T-dcpendcnt vari- 
ations of the total energies to be about 1.5 Ry/el./La, 
the Coulomb term in eq. [j] becomes very large, about 
14. However, even if the Coulomb energy represents an 
important correction in this case, there are several facts 
indicating that other energies are involved, and that the 
response of the system (during a self-consistent, spin- 
polarized calculation) will prevent the system to absorb 
all of the energy. First, an exchange splitting £ which 
is not uniform over all sites will modify the energy. The 
spin-polarized calculation shows that £ is finite (3-5 mRy) 
only close to La. Second, although the amount of charge 
transfer due to £ and fc^T is similar, there are long tails 
of the Fermi-Dirac distribution which make the radial 
dependencies of Aq^ and AQ(T) different and make Uq 
non-transferable from one case to the other. Third, the 
effective derivative of the La-DOS, N' La (e,T), (taken as 
f(-&)N' f (e)de), is T-dependent, because the bare DOS 
does not have the same large derivative over a wide en- 
ergy range. The mechanism of charge transfer and mag- 
netism is therefore expected to depend on T. This is con- 
firmed in the spinpolarized results. Usually, by raising 
fcflT, the magnetic moment goes down. But a few cases 
are found where the moment is largest for ksT about 3 
mRy. The reason is that the sharpest DOS structure is 
a bit beside the position of E^ in that case. 

Calculations of total energies converge slowly and a 
comparison between paramagnetic and spin-polarized re- 
sults are not made in all cases. From the converged cases 
it is seen that the total energies for magnetic solutions 
are lower than the non-magnetic ones. This difference 
is largest when the moments are largest, 10-14 mRy /cell 
for a total moment of 0.10-0.19 /j,b- For the cases with 
moments smaller than 0.05 (1b, the difference is only 0-2 
mRy. 

Two sets of spin-polarized calculations were done for 
the 270-site cells, both with 4 and 10 k-points. One used 
the same linearization energies, which for the 80-site cell 
reproduced well the result with the large basis. It is 
the paramagnetic DOS functions from these parametriza- 
tions which are shown in fig. ||. The DOS values at Ep 
and magnetic results are displayed in Table 3. The other 
set use different linearization energies, which for the para- 
magnetic 80- and 270-site calculations lead to increased 
DOS and derivative of the DOS. The spin-polarized re- 
sults using 10 k-points show a much larger moment for 
the latter case compared to the first set. The reason is 
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that there is a region of a rather fiat (or less steep) DOS 
close to Ej? in the first case, as can be seen in fig. [|. 
The bold line of the DOS from the 270-site cell is seen 
to follow the thin line of the 80-site cell, but only within 
a few mRy of E^?. Further below Ep the curve is again 
steep, as is the case for the other set of calculations for 
the 270-site cell. The result is that the charge transfer 
is lower in the former set of calculations. However, this 
difference between the two sets of calculations diminishes 
at larger T when the steep region of the DOS a bit aside 
of E^ becomes involved. This is consistent with the T- 
dependence of the moment, since the largest moment in 
the first set of calculations shown in Table 3, is when 
k B T is 4 mRy. At low T (k B T ~ 0.5 mRy), the moment 
decreases further. The T-dependence of the moment is 
in general opposite for the other cases, meaning that the 
moments become very small when k B T ~ 3 mRy. For 
the other set and in most cases with the 80-site calcula- 
tions the moments are largest for the smallest k B T ~ 0.5 
mRy. 

When the virtual crystal approximation is applied to 
this case, to remove 0.5 electrons from the system (and 
0.5 nuclear charges from the Sr sites), one finds that E^ 
falls closer to the steep left-hand side of the La-DOS, cf. 
fig |^. The objective with this model calculation is to 
see if the new conditions make the moment to increase, 
as was expected from the previous discussion. It should 
also approach the situation with two-component doping, 
where Al-doping (and a rigid-band shift of E^) in ad- 
dition to La-doping could influence the magnetic state. 
The result is that the spin-polarized results converge to- 
wards a larger moment, about 0.09 /i B , but now only 
for low T. At large k B T ~ 3 mRy, the moment shows a 
steady decrease. These T-dependent results corroborate 
the finding that there is a close connection between sharp 
DOS structures, charge transfer and magnetic moments. 

The results in Table 1 and 3 show that the moments 
per impurity are comparable for the small and the large 
cell, about 0.1 [i B per cell. The correlation between 
large derivative of the DOS and the moment, is seen in 
the results, both from the variation of the number of 
k-points and from the different choice of linearization. 
The largest moment, 0.25 /is /La is only slightly larger 
than the largest moment in the 80-site calculation, 0.19 
/is/La. This is despite the much larger derivative of the 
DOS from the results with 270 sites. As was discussed 
above, the reason could be that the large cell has a com- 
paratively larger undoped bulk, which opposes polariza- 
tion. 



C. Discussion. 

The results in ref.0 for an 80-site cell, showed a stable 
ferromagnetic state, and the moment remains centered 
around the La impurity. Together with the paramag- 
netic results, one can conclude that the magnetism is 



correlated with the derivative of the DOS on La, and 
that the magnetic state is partly stabilized by a charge 
localization on the La impurity. This correlation becomes 
more evident in the continued calculations using smaller 
basis, which are modified in order to create different con- 
ditions for large or small charge transfer effects. The 
results from the 270-site calculations show that the con- 
ditions for charge transfer to the La impurity are rein- 
forced, because the La-d band becomes even more narrow 
and localized due to the increased La-La distance in the 
large cell. The magnetic moments per La-impurity are 
comparable in the two cells though, probably because of 
the larger "bulk" between the impurities, which is un- 
doped and counteracts magnetism better in the large 
cell. On the other hand there is no magnetic state for 
the In or Al impurity, for the Sr-vacancy or for the un- 
doped SrBg. This is despite a certain pre-condition for a 
charge transfer from B towards remaining Sr in all these 
cases, but apparently it is too weak to produce an itiner- 
ant magnetic state, at least in the 80-site cell. The weak 
transfer towards Sr is present even in the La-doped cases, 
and therefore it reduces the (larger) transfer towards La. 
This competition could reduce the amplitude of the La- 
moment and even prevent magnetism if it became larger. 
One could speculate that Al substitutions on Sr sites play 
a role of a charge reservoir in this situation. Namely, as 
was noted earlier, an Al atom on a Sr-site makes no large 
effects on the DOS except for a rigid-band like down- 
ward shift of Ej? . No explicit calculations of having both 
La- and Al-substitutions were done, but taken separately 
these results show impurity band and rigid-band behav- 
ior, respectively. A varying Al content in a La-doped 
material should therefore modify the position of E^ rel- 
ative to the narrow La-d band, and create more or less 
optimal conditions for the charge transfer around the La- 
impurity. By tuning with Al-doping to find the largest 
DOS derivative on the La-band, one should be able to 
optimize the magnetic moment. 

The non-rigid band like electronic structure around the 
dilute La impurity shows up in the distribution of the 
ferromagnetic moment. According to these results, it is 
striking that La and the nearest sites takes the major 
part of the total moment. B-sites remain essentially non- 
polarized. The La-site acts as an impurity and will be a 
perturbation for excitons. Ferromagnetism has been pro- 
posed to be caused by the excitonic mechanismliU, since 
the undoped band structure show some favorable condi- 
tions for this. However, it is seen that with La-doping 
there are large modifications of the undoped bandstruc- 
ture, and a rigid-band assumption is not at all appropri- 
ate. Rigid-band like conditions occur for In-, A1-, to some 
extent for "vacancy" -doping, with Ef lower than for the 
undoped case. Theories of spontaneous ferromagnetism 
of the dilute electron gas have also been proposedE^Q. 
However, an unrealistically large electron-gas parameter 
r s is necessary for this (> 30), whereas the calculated r s 
is smaller than about 1.5. 

Some experiments can be suggested to verify if the La- 
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band is localized and if the moment is concentrated close 
to the La. Nuclear magnetic resonance experiments can 
detect a magnetic moment. But the magnetic moment 
is small and for a total moment of 0.19 \ib one finds 
a La contact field of less than about 40 kG. The spin- 
splitting is not exceeding 5 mRy within the La valence 
and core levels, which is small to be observed directly by 
spectroscopic methods. Furthermore, the low La concen- 
tration for the real magnetic system might be a problem 
for the intensity. The moment is mainly localized to the 
impurity, with weak interactions with the exterior. This 
situation resembles that of a Kondo impurity, but due to 
the small moment it is questionnable if this can produce 
the typical manifestations on specific heat and resistiv- 
ity at low T. Instead of looking for the moment directly, 
one can propose to use spectroscopic methods like pho- 
toemission, soft X-ray emission spectroscopy (XES) or 
soft X-ray absorption (SXA), to determine the degree of 
localization of the La-band. This can be done at larger 
La concentrations, before magnetism sets in. Such mea- 
surements could reveal if the DOS shows a localized La- 
band similar to the DOS calculated here for x ~ 0.125 
or even more at x ~ 0.037, or if the DOS is more like a 
rigid-band image of the undoped DOS. However, the low 
concentration of La can be a problem in some cases and 
lead to a inconclusive result. Soft x-ray photoemission 
(XPS), for instance, gives a often a good picture of the 
total DOS modulated by matrix elements, but if the La 
concentration is very small one can miss the interesting 
part regarding the localized impurity state. XES and 
SXA should be more useful, since they allow for a local 
projection of the bandstructure by the involvement of a 
core hole. A La core state can be used to probe the oc- 
cupied or unoccupied La valence states in XES or SXA, 
respectively. The signal might be weak because of few 
La sites, but the spectrum has the advantage of being 
a selective ^-projection of the local La-DOS. This means 
that the signal from La should show a narrow band, while 
a wide unstructured spectrum (similar to the total DOS 
for the undoped material) is expected if a core level from 
Sr or an In-impurity is used. Dipolar matrix elements do 
some reshaping of the pure DOS functions, but they can 
be estimated quite easily. The question about relaxation 
around the core hole is more complicated. According to 
the so-called final-state p4e, one observes the spectrum 
of the relaxed, final stateRil. If so, one should prefer XES, 
since the relaxed state is the unperturbed ground state, 
and the complication of a new calculation including a 
core-hole for SXA, can be avoided. 

Theoretical and experimental efforts have been done to 
determine the degree of baud, overlap at the X-point in 
the undoped hexaboridesEn3'EjE3. Calculations using gra- 
dient corrections to the LDA or using ihe GW method, 
find different values of the band overlapEm. If the overlap 
is absent and no doping, the material will be a small-gap 
semiconductor. The excitonic mechanism can be sensi- 
tive to this, and with pressure or measures for changing 
the z-parameter one can think of tuning the material for 



optimal properties. However, the mechanism for mag- 
netism that we have discussed here is not sensitive to the 
band structure of the undoped case. It is very sensitive 
to the degree of localization of the La impurity band. If 
a band calculation gives a too small localization, a wider 
band with small derivatives of the DOS, then the moment 
goes to zero. So, unless there is an unknown connection 
between the band gap for the undoped band structure 
and the localization of the La impurity band, we believe 
that the exact shape of the bands in the undoped system 
is irrelevant for the mechanism for magnetism around a 
La- impurity. But, since we have seen that the results 
are sensitive to features of the impurity band, it also is 
probable that lattice relaxation around the impurity, in- 
cluding variations of the z-parameter, will have an effect 
on the magnetic moment. 

This work show that the effect of dilute La-doping on 
the band structure of SrBg can lead to the appearance of 
a magnetic moment around the impurity. Unfortunately, 
temperature variations cannot be studied with the same 
level of sophistication. The experimental observation of 
magnetism is interesting, but also its temperature depen- 
dence is unususal as magnetic ordering can persist up to 
about 600 Kb. The Fermi-Dirac smearing is only one part 
of the T-depencencies. It is known that thermal struc- 
tural (vibrational) disorder is important already at 200- 
300 K. Delicate changes of band structure and the fine de- 
tails of the DOS due to thermal disorder can be essential 
for finding a correct description of the physical properties 
of a materialE3. This, disorder in the distribution of im- 
purities, disorder of the local moments, non-collinearity 
and orbital magnetism are not considered here, only the 
electronic excitations given by the Fermi-Dirac distribu- 
tion. In general, the calculated magnetic moment dis- 
appears for UbT in the range 2-4 mRy (300-600 K), but 
not always. One case with the large cell had the largest 
moments for quite large T, which could be related to 
a sharp DOS structures a bit beside the optimal posi- 
tion for Fip. If this mean-field picture is sufficient, one 
would then expect a wide range of ordering temperatures 
among different samples, like different samples showing 
different amplitudes of the moments. The effects of spa- 
tial irregularities of the impurity distribution of La or 
other dopants have not been calculated explicitly. But 
from the separate results with different relative positions 
between and the localized La band, one can imagine 
that some samples could be apt for large moments at low 
T, while others with lower moments could resist better 
at large T. 

Finally, the question is whether this mechanism for 
magnetism is rare and exotic, or if systems other than 
the hexaborides can show similar effects. From eq. |l| one 
should have a small total DOS, but large DOS derivative 
of a sub-band. Pure compounds with wide, hybridized 
bands can hardly be interesting cases. Dilute doping in 
a low DOS material seems to give optimal conditions, 
because N v /N — > 1, and an impurity band can localize 
and become narrow for large distances between impuri- 
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ties. Not all combinations between impurities and host 
material will do, because the hybridization can make the 
impurity band to merge into the host bands, as for In in 
the hexaboride. The requirement of a large Uo with the 
correct sign, will also reduce the number of candidates 
for magnetism of this kind. Preliminary calculations for 
metal impurities in some oxides show similar results as 
for the La doping here, although there is-a, sizable polar- 
ization in the proximity of the impurityEij. 



IV. CONCLUSION. 



band" dopants like Al. It is probable that other systems, 
with localized bands from selected dilute doping in low 
DOS materials, can be found to show similar type of weak 
magnetism. 
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We have shown that ferromagnetism occurs in the su- 
percells of LaSr7B4§ and LaSr26Bi62- By doing the cal- 
culations with different number of k-points and different 
basis functions, it has been possible to show that the 
actual derivative of the local DOS at Ep, which is neces- 
sary for charge transfers, scales with the stability of the 
magnetic moment. This demonstrates that the mecha- 
nism for ferromagnetism is not exclusively found in the 
exchange energy. The necessary condition of a localized, 
non-constant La-DOS near Ep, appears to be enforced as 
the doping decreases, and the mechanism of charge trans- 
fers persists at lower doping levels. Thus, it is likely that 
the observed W£ak ferromagnetism for dilute La-doping 
in hexaboridescm, is a band effect of the kind that we 
have described. This conclusion is based on results from 
first-principle LMTO calculations. The moments are of 
the order 0.05-0.1 /is/La for the calculations with the 
"best" choice of basis function and the largest number 
of k-points. However, the fact that the results can be 
"boosted" by choosing somewhat extreme linearization 
energies, brings up the problem of describing a very sen- 
sitive band structure with a linearized method. Most 
band results (total energy, Fermi surface, and the general 
shape of bands and DOS functions) usually show only 
a moderate sensitivity to the choice of basis or k-point 
convergence. But here, although the visual inspection of 
the DOS show no large differences, there are small differ- 
ences in the localization of the impurity bands depending 
on the choice of basis, which turn out to be important for 
the size of the magnetic moment. Therefore, one must 
emphasize that the calculations show the likelihood of a 
ferromagnetic state, but that one can not be affirmative 
for the absolute values. This situation may reflect the 
reality, since different samples are reported to have dif- 
ferent moments (at most 0.32 ^s/La) , -while some show 
only large paramagnetic susceptibilitiesB. Real materials 
have some randomness in the distribution of the impuri- 
ties, while the calculations are made for ordered super- 
cells. Clustering and other lattice imperfections can be 
important for the properties if they lead to a rigid-band 
like modification of the band structure, or if the impurity 
states become more localized. From the balance between 
such effects, it is suggested that optimal conditions for 
magnetism can be found indirectly by addition of " rigid- 
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FIG. 2. Paramagnetic, total and partial La-DOS (includ- 
ing the DOS from the 6 closest empty spheres) near Ef for 
LaSr7B48 and LaSr2eBi62 as obtained from tetrahedron inte- 
gration for 20 and 10 k-points, respectively. The zero energy 
is at Ef in both cases. Note that the DOS is per formula unit 
(SrBe), so that the total DOS is divided by 8 in the former, 
and by 27 in the latter case, respectively. Thus, the partial 
La DOS is much larger in the large cell (about 55 per La) 
compared to in the small cell (about 20 per La). 
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FIG. 1. Paramagnetic, total DOS for LaSr7B48 (bold 
line), InSr7B4g (broken line) and undoped SrgB4g as obtained 
from tetrahedron integration for 20 k-points using the small 
basis. The zero energy is at Ef for each case. Note that 
In-doping leads to a downward rigid-band like shift of Ef on 
a DOS which resembles much that of the undoped material. 




TABLE I. Paramagnetic DOS values (states/cell Ry, and 
states/cell Ry ), and magnetic moment on La and total (/j,b), 
as function of the number of k-points for the 80-site unit cells. 
To the impurity site is counted also the contributions from the 
6 closest empty spheres. Values within parentheses are from 
the calculations with the reduced basis. 



k-points 


N(E F ) 


N La 


N' La 


m La 


mtot 


10 


52(54) 


23(24) 


830(830) 


0.10(0.07) 


0.15(0.12) 


20 


48(48) 


19(21) 


730(840) 


0.07(0.06) 


0.10(0.09) 



TABLE II. Differences in valence charge (el. per site) 
between doped and undoped SrB6-supercell for selected sites, 
calculated from 20 k-points in the small unit cell. The charge 
within the 6 closest empty spheres is assigned to the impurity 
site. There is not much variation among non-equivalent B 
sites. 



Basis 


Impurity-site (La or 


Sr) Sr near La Sr far from La B( 


Large 


1.19 


0.005 0.010 


Small 


1.15 


0.006 0.009 



-0.005 
■0.001 



TABLE III. Paramagnetic DOS (states/cell Ry), and 
magnetic moment on La and total (/is), as function of the 
number of k-points for the 270-site unit cells. To the impu- 
rity site is counted also the contributions from the 6 closest 
empty spheres. The derivatives of the La-DOS are roughly 4 
times larger than in the results from the 80-site cell, but the 
non-linearity is large and no precise numbers can be given. 
The values within parentheses are for a different linearization, 
as is explained in the text. The magnetic moments are cal- 
culated at low fesT ~0.5 mRy, except the one for 10 k-points 
without parenthese, which is for fcsT=4 mRy, see the text. 



k-points 


N(E F ) 


N La 


m La 


mtot 


4 


230(220) 


90(80) 


0.13(0.17) 


0.19(0.25) 


10 


165(270) 


55(85) 


0.04(0.09) 


0.05(0.17) 



-0.05 0.05 

(E-E F ) (Ry) 
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